i Boo (0.2678) 




HT Vir (0.4077) 




VW Pic (0.4258) 



Excl 
_i I 



7.05 
_i I 



V445 Cep (0.4488) 



6.87 
_i I i_ 



V867 Ara (0.4937) 



Excl 
_i I 



7.43 
_i I i_ 



V918 Her (0.5748) 



T 1 1- 



7.40 
_i I i_ 



RR Cen (0.6057) 

















V 




7.40 



LL Vel (0.6192) 

I I 



6.81 
_i I i_ 



VW Cep (0.2783) 



V759 Cen (0.3940) 



\ «■ * f" \ «■ 



y ■ 



7.46 




KR Com (0.4080) 



V566 Oph (0.4096) 




AW UMa (0.4387) 



_i I I i_ 



6.92 
I i_ 



CN Hyi (0.4561) 




KS Peg (0.5021) 





^'--^ ^- 






5.46 


V353 Peg (0 


5846) 








7.49 



IS CMa (0.6170) 




V535 Ara (0.6293) 
-I 1 1 1 1 1— 



_i I I i_ 



7.22 
I I 



V972 Her (0.4431) 



6.72 
_i I I i_ 



MW Vir (0.4931) 





1./. - • 


• " -w.«.w - 




7.01 



AA Cet (0.5362) 




eps CrA (0 


5914) 














4.84 



V851 Ara (0.6174) 
-I 1 1 1 1 1— 



_l_ 



6.93 
_i I I i_ 



S Ant (0.6483) 
-I 1 1 1 1 1- 




VIOIO Oph (0.6614) 




KP Peg (0.7272) 



7.10 
_i I i_ 



VI 130 Tau (0.7989) 



Excl 
_i I 



6.64 
_i I i_ 



TV Pic (0.8520) 




DW Boo (0.9244) 









7.40 



V449 Aur (0.7036) 



V2082 Cyg (0.7141) 




PP Hya (0.7523) 



FO Vir (0.7756) 




V2388 Oph (0.8023) 



V335 Peg (0.8107) 




ES Lib (0.8830) 



IW Per (0.9172) 




TU Hor (0.9360) 



GK Cep (0.9361) 




HL Dra (0.9443) 

I I I I I I 

** * ^ J 1 ^ s 



Excl 7.34 
_i I I I I i_ 



DX Aqr (0.9450) 



T 1 1 1 1 1- 



Excl 
_i I i_ 



6.45 
I i_ 



BB Sol (0.4765) 
-I 1- 



ER Vul (0.6981) 

Tl 1" 



0.10 



0.08 - 



0.06 - 



0.04 - 



0.02 - 



0.00 



0.02 



DX 



1010 



ES 



AA 



#F0 



KP 



1 



_L 



_L 



I A I 



_L 



0.00 -0.05 -0.10 -0.15 -0.^ 

3.0 



0.10 



0.08 - 



0.06 - 



0.04 - 



0.02 - 



0.00 
0.02 



DX 



1010 



ES 




TV 



TU 



-0.020 -0.010 0.000 0.010 0.020 



U.U 




1 1 1 

+ + - 

+ ^^^^^ • +^ 


0.2 




/ ^ + + o 


> 0.4 






1 

O.b 




• 


0.8 


/• 
/ • 




1.0 




1 1 1 



0.2 0.4 0.6 0.8 1.0 

Period (d) 



4 



0.6 



- 1- - BW. 



- 0.4 



o 

CO 

I 

O 

a. 2 

in 
I 

O 



£ 1 



I 
I 

T 
I 

1 



5x 



- 0.2 



0.0 







0.2 0.4 



1"": 
- -I — : - 

- ■ ^ ■ I - [■■■-^■■■| - T 



0.6 0.8 

P(d) 



1.0 



The 7.5 Magnitude Limit Sample of Bright Short-Period Binary Stars. I. 
How Many Contact Binaries Are There? ^ 



Slavek M. Rucinski 

David Dunlap Observatory, University of Toronto 
P.O.Box 360, Richmond Hill, Ontario, Canada L4C 4^6 

rucinskiOastro . utoronto . ca 
ABSTRACT 

A sample of bright contact binary stars (W UMa-type or EW, and related: with 
j3 Lyr light curves, EB, and ellipsoidal, ELL - in effect, all but the detached, EA), to the 
limit of V'^^^ = 7.5 magnitude is deemed to include all discoverable short-period (P < 1 
days) binaries with photometric variation larger than about 0.05 magnitude. Of the 32 
systems in the final sample, 11 systems have been discovered by the Hipparcos satellite. 
The combined spatial density is evaluated at (1.02 it 0.24) x 10^^ pc^^. The Relative 
Frequency of Occurrence (RFO), defined in relation to the Main Sequence stars, depends 
on the luminosity. An assumption of RFO ~ 1/500 for My > +1.5 is consistent with the 
data, although the number statistics is poor with the resulting uncertainty in the spatial 
density and the RFO by a factor of about two. The RFO rapidly decreases for brighter 
binaries to a level of 1/5,000 for My < +1.5 and to 1/30,000 for My < +0.5. The high 
RFO of 1/130, previously determined from the deep OGLE-I sample of Disk Population 
W UMa-type systems towards Baade's Window, is inconsistent with and unconfirmed 
by the new results. Possible reasons for the large discrepancy are discussed. They 
include several observational effects, but also a possibility of a genuine increase in the 
contact-binary density in the central parts of the Galaxy. 

Subject headings: stars: close binaries - stars: eclipsing binaries - stars: variable stars 



1. THE CHANGING VIEWS ON THE SPATIAL DENSITY OF THE W UMA 

SYSTEMS 

Our views about the spatial density of contact, W UMa-type binaries have meandered con- 
siderably during the recent half century. Shapley (1948) was the first to point out that "W Ursae 



^Based on data from the Hipparcos satellite mission and from the David Dunlap Observatory, University of 
Toronto. 
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Majoris variables ... are not only the most numerous of eclipsing stars ... but ... more numerous 

than all other variable stars combined". From the numbers given by Shapley, one can infer that 
he estimated that perhaps as many as about one percent of all solar-type stars are W UMa-type 
variables. 

Several attempts have been made to put the statement of Shapley in a quantitative way, with 
diverse results. Popov (1964) estimated the local density at 2 x 10^^ pc^^, while Kraft (1967) 
estimated it to be some 20 times lower, 10""^ pc~^. Then van't Veer (1975) found a much larger 
number, 11 x 10~^ pc~^, which corresponds to about one percent of all stars being W UMa-type 
binaries. Subsequently Duerbeck (1984) derived a value of 2± 10~^ pc~^, equivalent to the Relative 
Frequency of Occurrence (RFO) of ~ 0.001 (or 0.1 percent), and this value remained in popular 
use for some time. A similar, simplified, and thus rather convincing (based on naked eye objects) 
estimate of the frequency was given by Rucinski (1993), RFO ~ 0.0005 — 0.002, that is one W UMa 
binary per 500 - 2000 ordinary stars. 

The newest investigations, based on much larger statistical samples, suggested that the RFO 
is perhaps as high as suggested by Shapley, i.e. some five to ten times higher than estimated by 
Duerbeck (1984). During the last two decades, several EW systems have been discovered in old 
open clusters (Kaluzny &; Rucinski 1993; Rucinski &; Kaluzny 1994; Rucinski 1998b). A combined 
approximate estimate of Rucinski (1994) for a few old open clusters gave the RFO of ~ 275 ± 75 
ordinary dwarfs per one EW system. Later, in a combined analysis of several open clusters, Rucinski 
(1998b) showed that the RFO evolves with the age of the stellar system; the number of EW systems 
increases from the level of one per a thousand at the age of about 0.8 Gyr to the level of one per 
some 200 - 300 dwarfs at the age of 5 - 7 Gyr. The RFO was found also very high in globular 
clusters Rucinski (2000); however. Population II contact binaries are of no relevance to the present 
study which concentrates on the Population I objects, primarily in the solar neighborhood. 

The highest RFO was estimated for the Galactic Disk, as probed by the narrow, deep pencil 
beam of the OGLE-I survey in the direction of Baadc's Window (Rucinski 1995, 1997a,b, 1998a,b). 
This finding was in a basic accordance with the view of the increased RFO from the age of the 
oldest open clusters to the age of the Galactic Disk Field of about 11 Gyr (Binney et al. 2000). 
The great advantage of the OGLE-I result over all previous spatial density estimates was not 
only in the size and uniformity of the sample (a few hundred contact systems discovered in the 
same survey), but - mainly - in that the RFO estimate was based on the volume-limited samples, 
complete to My ~ 5.5 to the distance of 3 kpc, and (with better statistics for brighter systems) to 
My — 4.2 to the distance of 5 kpc. Large numbers of binaries in the OGLE-I sample permitted 
to evaluate the RFO on the per-My-bin way, in place of the previous estimates based on data 
averaged for all accessible spectral types, between late A to early K. Such comparisons through 
the luminosity function led to a statistically well established - and high - apparent frequency at 
the level of ~ 1/130 (Rucinski 1998b). A clear drop in numbers for brighter, longer-period contact 
systems, with a sharp cut-off at P ~ 1.3 — 1.5 days, was also noted (Rucinski 1998a,b). 
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The OGLE-I result would suggest a return to Shapley's value of the RFO at the level of one 
percent or higher. However, this would mean a rather large disagreement with the relatively firm 
result of Duerbcck (1984) for the solar neighborhood. Can we exclude a systematic error in the 
OGLE-I result? Even if based on the best number statistics, the OGLE-I may have systematic 
biases. After finding that the new 7.5 magnitude-limit sample indeed does not agree with the 
OGLE-I result, we attempt to give possible reasons for the discrepancy in Section 11. 

Sections 2-7 describe construction and properties of the sample used in a new estimate of 
the spatial density. This estimate is through the luminosity function, as described in Section 8. 
The related "period function" (the number of systems per unit of volume, per period interval) 
is described in Section 9. The amplitude distribution for the sample is discussed in Section 10. 
Section 11 discusses the major discrepancy in the density spatial estimates between the current 
sample and the OGLE-I sample in the direction of Baade's Window (Rucinski 1998b). Section 12 
looks into the future of the spatial density estimates for contact binaries. Section 13 summarizes 
the results of the paper. Appendix A gives brief descriptions of the individual binaries, including 
those which have been excluded from the sample in the last stages of its definition. 

2. CRITERIA FOR SELECTION OF THE 7.5 MAGNITUDE SAMPLE 

A well defined volume-limited sample would be the proper way of estimating the spatial density 
of the EW systems. The Hipparcos sample immediately comes to mind because - in principle - 
the parallax data can be directly used to place systems in space and thus lead to straightforward 
estimates of the numbers per unit of volume. The Hipparcos sample of bright stars is also currently 
the most complete as far as stellar variability is concerned, down to the amplitude levels perhaps 
0.02 - 0.03 magnitude, and certainly to amplitudes larger than 0.05 magnitude. As we will see, 
in the final sample that we considered of 32 close binary systems, 11 systems are new Hipparcos 
discoveries, so that the increase in numbers is significant and must be considered in the spatial 
density estimates. However, the completeness of the Hipparcos sample is in fact limited to the bright 
stars, to y ~ 7.5. The much deeper Hipparcos sample of the W UMa-type systems discussed by 
Rucinski &; Duerbeck (1997) was not complete in any sense as it relied on the input lists of all then 
known systems which was provided to the Hipparcos project when the mission was programmed. 
Thus the list of Rucinski &; Duerbeck cannot be used to evaluate the space density of the W UMa- 
type systems. 

Of the two advantages of the Hipparcos data for bright stars, the availability of parallaxes 
and the completeness for photometric variability detection, the latter appears - paradoxically - to 
be more important. As is well known, Hipparcos obtained the first set of reliable trigonometric 
parallaxes. But they still had relatively large random errors, as can be illustrated by the statistics 
of parallaxes for our final sample. While the maximum value of the parallax (44 Boo) is 78.39 mas 
(milliarcsec) , the median value for the sample is only W = 12.40 mas, with the smallest value being 
4.37 mas. This should be compared with the median of the mean standard errors a{Tr) = 0.85, with 
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the smallest and largest errors being 0.41 and 2.36 mas^. Recalling that, as was shown by Lutz 
& Kelker (1973), large systematic errors in distances start occurring for parallaxes smaller than 
TT ~ 7 a{Tr), we cannot use the parallax data to select the sample. 

Thus, knowing that construction of a volume-limited sample still remains impossible, we try 
to construct a complete magnitude-limited sample and then use it to obtain the spatial density 
estimates in intervals of My. With this goal in mind, we constructed the 7.5 magnitude-limit 
sample, in the following way: We selected the short-period binaries by merging data from the two 
basic sources, the Hipparcos Catalog (ESA 1997, HIP) and the General Catalogue of Variable Stars 
(Kholopov et al. 1999, GCVS). The selected binary stars had the orbital period shorter than one 
day and were brighter than 7.5 magnitude at maximum light. Since the number of systems rapidly 
increases with the limiting magnitude, the magnitude cutoff was particularly carefully observed, as 
described in Section 3. The version of the GCVS which was used is available only by the Internet^ 
and consists of Volumes I - III of the Fourth Edition of 1988 merged with the Name Lists Nos. 67 
- 74. We added to this the Name List 75 (Kazarovets et al. 2000). This material was considered 
complete in the Fall of 2001. To limit the sample to the nearby stars, we removed stars with the 
HIP parallaxes smaller than 2.5 milli-arcsec (mas), i.e. more distant than 400 pc, and, thus - with 
our magnitude limit - brighter in absolute sense than My = —0.5. 

The sample, at this stage consisting of 113 variable stars, underwent further scrutiny for 
presence of pulsating stars, spurious variables or variables with no orbital periods (formally zero 
days), but with single minima (probably detached eclipsing binaries) observed by Hipparcos. All 
relevant literature information was used, particularly on S Set and 7 Dor pulsating stars, leading 
to a reduction in the sample to 46 stars. Of particular interest are the following items: 

• Three variables are included in the HIP Catalog as EB systems, but appear to be pulsating 
variables. This result is based on our own spectroscopic data obtained at the David Dunlap 
Observatory (Rucinski et al., private communication). These are the small amplitude (0.04 
- 0.05 mag.) variables: CC Lyn (HD 60335, HIP 36965, V"*"^ = 6.40), FH Cam (HD 63383, 
HIP 38900, = 6.88) and CU CVn (HD 120349, HIP 67357, F"^"^ = 7.53). The two 

former show sharp-line spectra with radial-velocity changes within 10 to 20 km s~^ and 8 to 
14 km s~^, while the third (whose radial velocity variations remain to be analyzed) shows 
rotationally broadened spectra with V sini = 155 it 10 km s^^. In all three cases the periods 
appear to be half of the HIP periods, i.e. 0.177311, 0.13624 and 0.067834 days. Note that 
without spectroscopic data these three stars would be classified and included by us among 
the ellipsoidal binaries (ELL). The final sample of ELL variables consists of 14 objects so that 



^Four binaries with the parallax errors larger than 2 mas, HT Vir, V867 Ara, DX Dra and AA Get, arc rricrnbors 
of triple systems. The three former have been eliminated from the sample for other reasons so that only AA Get 
remains in the sample. More about AA Get in Section 7, where the very low value of the parallax and its large error 
are specifically addressed. 

^ http : // www. sai .msu.su /groups /cluster /gcvs / gcvs /iii/ 
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removal of the three stars was significant; these numbers can be also taken as an estimate of 
the possible level of the contamination of the sample by the pulsating variables. 

• While we excluded all detached binaries from our sample, we scrutinized all which fell within 
our period and brightness limits in order not to miss relevant stars which could be incorrectly 
classified. Among the five detached binaries in the HIP catalog with EA designation, fulfilling 
our brightness and period criteria, three binaries: MX Hya (0.7960 d), UU Psc (0.8417 d) and 
V772 Her (0.8795 d), do not have reliable HIP light curves, and only two have reasonable HIP 
light curves: BB Scl (0.4765 d) and ER Vul (0.6981 d). We analyzed the latter together with 
the contact binaries, mostly as a check of our procedures. It should be noted that ER Vul 
frequently appears in various catalogues under an incorrect entry of EW (we expand on the 
naming conventions EW, EB, ELL and EA in Section 5). 

• We excluded two, most-likely EA binaries from consideration: BR Ind and V2300 Oph. While 

BR Ind appears in the HIP Catalog with the period of 0.8928 day, its period is most likely 
twice that long. The same possibility applies to the extremely small amplitude variable 
V2300 Oph (Jerzykiewicz 1993) with the nominal period of 0.9071 day. 



3. Hp AND V MAGNITUDES AND THE 7.5 MAGNITUDE LIMIT 

Our selection of the sample and assignment of variability types was based on the Hipparcos 
light curves in the instrumental Hp system. Since the number of objects quickly increases with the 
limiting magnitude, it is important to note that the cutoff at the maximum-light 7.5 magnitude is in 
the standard V magnitude bandpass. For each Hipparcos light curve, the maximum Hp magnitude 
was determined. The values are given in in individual panels of Figures 1 and 2. The order in these 
figures is the same as in Table 1 and follows the increasing orbital period. 

The values of i/p have been transformed into V"^""^ using the table of the differences V — Hp 
as a function V — I, as given in Table 1.3.5. in the HIP Catalog Explanatory Books. The V — I 
indices were in turn estimated from the mean B — V values, as derived from the Tycho Bt and 
Vt data in the Tycho-2 catalogue (Hog 2000, TYC2), using the appropriate linear transformations. 
Since most of the objects do not have ground-based photometric indices, this source of the average 
B — V data assures some uniformity for inter-comparison. However, the resulting values of V"^"'^ 
may be uncertain because of the two-step interpolation of V" — Hp. For that reason, whenever 
available, ground-based values of y™-"^ have been used in preference to the values determined from 
jjmax_ rp^^ systems, V759 Cen and V566 Oph, have i^^"^ below the cutoff limit, but their 
transformed place them at V'^"'^ < 7.5, confirming the available literature data. 

Table 1 contains the relevant data for the systems of the sample. Notes on individual objects 
are given in the Appendix A. Of the 43 objects listed in the table, two (BB Scl and ER Vul, listed 
at the end of the table) were known from the beginning to be EA systems. Six EA systems were 
later identified among binaries originally classified as EW and EB (the variability-type classification 
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is discussed in full in Sections 5 and 6) . Finally, three W UMa-type systems in triple systems have 
been rejected because only the combined brightness with their close companions was above the 
adopted magnitude limit; this is discussed in the next Section 4. 

The final sample of contact and related stars (EW + EB + ELL) consists of 32 objects. The 
sample does not contain many popular systems repeatedly observed over the years, among them the 
prototype of the class, W UMa. On the other hand. Table 1 contains 17 new Hipparcos discoveries, 
of which 11 have survived the selection scrutiny to be included in the final sample. These numbers 
best illustrate that a re-discussion of the current statistics, even for a magnitude-limited sample 
with a bright limit, was very much needed. 

4. THE THIRD LIGHT: CLOSE COMPANIONS 

Several among the binaries of the sample have speckle interferometry, close visual or spectro- 
scopic companions. The light contribution of a third component will be characterized here by the 
fractional light factor, /3 = L^/L"lf'^ , the ratio of the observed luminosity of the third component 
to the observed luminosity of the close binary at its maximum brightness. The values of (3 for 
individual cases were found from the literature, as given in the Appendix, usually after consulta- 
tion of the very useful studies of Tokovinin (1997) and Fabricius & Makarov (2000). We list the 
adopted values of ft in Table 1. We also give there comments on the angular separation (in arcsec). 
The final, corrected values of the maximum magnitude for the binary systems V"^"'^ have been 
computed using: V"^"'^{corr) = V"^"-^{obs) + 2.51og(l -|- (3). 

The current sample has been reduced in size by taking into account the presence of close "third 
components" in systems which - without allowance for the companions - were brighter than the 7.5 
magnitude, but which became fainter than this limit when the "third light" was subtracted. The 
three systems, included in Table 1, but excluded from any further statistical consideration were 
HT Vir, KR Com and V867 Ara. The downward shifts from the observed values of V"^"^ to the 
corrected ones were moderately small (0.41, 0.16, 0.24 mag.), but significantly large to cross the 
7.5 magnitude limit. We are particularly sorry to part with the system HT Vir which is one of the 
most interesting in the contact-binary world: It shows a large amplitude of light variations (when 
corrected for the third light), in perfect agreement with its large mass ratio, qgp = 0.812 it 0.008 
(Lu et al. 2001). 

Presence of a third light in some system leads to a "dilution" of the variability signal in the 
light curve with the corresponding decrease in size of the Fourier coefficients used for variability 
classification (see Section 6). To account for this effect, the light curve decomposition coefficients, 
expressed in light units, have been corrected for the third-light using: ai{corr) = ai{obs){l+l3). The 
light variation amplitudes (amp), measured in magnitudes, also change according to: amp{corr) = 
-2.51og[(l + (3) dex(-0.4 amp{obs)) - /?]. 
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5. VARIABILITY TYPES 

The goal of this paper is to evaluate the spatial density of contact binary stars. Three classes 
of binaries must be considered in this context, the EW or W UMa-type, the EB or Beta Lyrae-type 
binaries and the ELL or Ellipsoidal variables. The distinction between them is sometimes difHcult, 
especially when color curves are not available. 

The EW or classical contact or W UMa-type binaries are characterized by continuous light 
variations and equally deep eclipses. The latter property results from equal temperatures of compo- 
nents, in spite of their differing masses. The equality of the temperature and thus very good thermal 
contact can be taken as the best defining characteristics of the contact-binary class. Normally, EW 
binaries show a slight reddening at both eclipses due to the combined effects of the gravity and 
limb-darkening which dominate over a small temperature difference between the components. 

Obscrvationally, the class EB is defined as basically a modification of the EW class, but with 
unequally deep eclipses, leading to light curves of the P Lyrae-type. Since the differing depth 
of the minima is a component temperature-difference effect, the color curves show reddening and 
blueing in primary and secondary minima, respectively. Spectroscopists tend to call EB systems 
those that show lines of only one component. Physically, what is observed as the EB binaries is, 
unfortunately, a mixture of distinct classes of (1) contact binaries in poor thermal contact, (2) 
pre-contact or broken-contact semi-detached binaries and (3) semi-detached Algols, after the mass- 
exchange and mass-ratio reversal. Thus, the classification of EB may be occasionally misleading. 
The real, massive and evolved, /? Lyrae-type binaries (sometimes also called W Ser binaries), with 
the inverted mass-ratio and on-going accretion into a large accretion disk, do not exist among very 
short-period stars (there is no space for the disk). There are also indications that for periods shorter 
than one day, Algols arc extremely rare (Rucinski & Lu 2000). The most frequently occurring 
short-period EB objects appear to be binaries either just before establishing contact or in one of 
the broken-contact semi-detached stage, in any case, prior to the mass-ratio reversal, with the more 
massive component close to or at its critical Roche surface. 

The distinction between the EW and EB class in individual cases is frequently difficult even 
when radial velocity data for both components are available. In our radial-velocity program con- 
ducted at the David Dunlap Observatory, we have observed systems which have been originally 

called EB on the basis that only one set of spectral lines had been observed in previous programs, 
but for which we see now spectra of both stars, sometimes suggesting a different classification. 
Frequently, no directly available data can tell us if the star is an EW system with a broken thermal 
contact or an EB system with the more massive star filling the Roche lobe and an undersized sec- 
ondary. In both cases, the more massive component is the hotter one. Both types are evolutionary 
related with one preceding or leading to the other. Thus, we decided to use both classes, EW 
and EB, in a combined spatial density estimate. We note that Stepien et al. (2001) (also Stepien 
2001, private communication) suggested that contact binaries are basically semi-detached objects 
with the smaller component engulfed by the matter of the more massive, larger component. This 
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presumably dynamically-stable modification of the Shu model (Shu &: Lubow 1981) still requires an 
actual modeling effort. If this model would work, the distinction between the EW and EB classes 
would reduce only to the efficiency of the energy transport. 

The third class are the ellipsoidal, or ELL, variables. This class has no real physical meaning, 
as these are - in majority - EW and EB binaries seen at small orbital inclination angles. We 
simply have no way of assigning these low-amplitude variables to the two classes. Predictions on 
the numbers of low amplitude contact systems (Rucinski 2001) suggest that there should be many 
such variables and indeed the Hipparcos mission discovered majority of them among the bright 
stars. As discussed in Rucinski (2001), some ELL variables may be very low mass-ratio objects 
seen at large inclinations, but then total eclipses can be observed even for relatively low inclinations. 

Both catalogs used in the selection of the current sample, the HIP Catalog and the Variable 
Star Catalog, give the variability types. However, we found that these are frequently unreliable. 
This was particularly so for the large number of small-amplitude systems called EB in the HIP 
Catalog, most probably by some sort of an automatic classifier. Of course, these could be in fact 
EW or EB binaries, but the tendency to call all variables EB seemed to increase for very small 
amplitudes, exactly in the regime where usually nothing can be said solely on the basis of the light 
curve. For that reason, we used the name of ELL variables to signify all small amplitude systems 
with amplitudes smaller than 0.1 magnitude (more strictly, when \a2\ < 0.05). 

Our goal was to determine the density of contact binaries which we here identify with the 
combined EW + EB + ELL classes. The detached (EA) binaries have been eliminated early from 
consideration, mostly because the statistics of short-period detached binaries is practically impos- 
sible to establish at this time on the basis of the photometric data alone. Also, while reasonably 
secure predictions on the number of systems missed due to low inclinations can be made for stars 
filling the Roche lobes, mostly because of the relatively simple geometry (Rucinski 2001), such 
predictions are entirely impossible for detached binaries for which the parameter space is too large. 
In the final sample, we identified and eliminated the EA binaries through the Fourier analysis of 
the Hipparcos light curves, as described in Section 6. 

6. FOURIER ANALISIS OF THE LIGHT CURVES AND VARIABILITY 

CLASSIFICATION 

The light curves obtained by the Hipparcos mission in the Hp wide-band photometry were 
retrieved from the HIP database for all 43 close binaries with periods shorter than one day and the 
maximum brightness above the V = 7.5 magnitude limit. We show the light curves in Figures 1 
and 2. 

The classification and selection of the objects has been done by performing the Fourier analysis 
of the Hipparcos light curves in the same way as described originally in Rucinski (1997a,b). The 
Fourier coefficients are listed in Table 1. They have been determined by least-squares fits to 
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light curves expressed in light units (using the assumed H'p"'^ as the reference level), with errors 
estimated through a "bootstrap" re-sampling process. The median value of the mean standard 
error for all coefficients was 0.0018 (this value is identical for all coefficients for a given light curve 
because of the orthogonality), but with some errors - for poorly-covered light curves - reaching 
much larger values than the sample median. Note in particular that oi, which by definition should 
be negative, is positive for some systems. The positive values of a\ are invariably associated with 
large errors due to poor light curves or incomplete phase coverage, e.g. for V759 Cen and DX Aqr 
the mean standard errors of ai are 0.0045 and 0.0055. For some stars, the coefficients and their 
errors have been magnified when allowance for the third light was made (Section 4). Such corrected 
coefficients are shown in the figures described below; they can be eaily identified by the unusually 
large error bars in the figures. 

The separation of the EW and EB systems from the detached binaries, EA, was done in the 
02 — 04 plane. The "contact" line in Figure 3 cleanly separates the detached binaries from the 
EW/EB systems which are located below this line. The ELL systems in this plane are these 
binaries which arc not obviously EA and show small 02 (in the absolute sense), —0.05 < 02 < 0; 
this approximately corresponds to the variability amplitude < 0.1 magnitude. Six EA binaries have 
been identified and removed from the sample at this stage: VW Pic, V449 Aur, VI 130 Tau and 
GK Cep appear to consist of strongly tidally-distorted, detached components of similar brightness, 
while HL Dra and DX Aqr, appear to be detached, but with components of different temperatures 
(see below). 

Figure 4 shows the 02 — oi coefficient plane which can be used to separate the EW and 
EB systems. The first cosine term, oi, is sensitive to the temperature difference between the 
components. We used the criterion of a\ < —0.02, the same as in Rucinski (1997b), to assign 
systems to the EB class. There are only five EB systems in the sample: KP Peg, FO Vir, ES Lib, 
AA Get and VIOIO Oph. The other systems above the ai = —0.02 line, DX Aqr and HL Dra, are 
the two EA binaries with components of different temperatures. 

The third of the Fourier-decomposition diagrams (Figure 5) addresses the light-curve asym- 
metry, measured by the first sine term, 61. It is clear that the O'Gonnell effect (Davidge &l Milone 
1984), of the first light-curve maximum being higher, is clearly present in VlOlO Oph and ES Lib. 
The type of the correlation between oi and hi was quite prominently visible in the OGLE-I data 
for EB binaries (Rucinski 1997b). It was interpreted as a result of accretion phenomena in a 
semi-detached configuration with the more massive star filling the Roche. 

7. THE PERIOD-COLOR RELATION AND THE ABSOLUTE MAGNITUDE 

CALIBRATION 

The period-color relation is shown in Figure 6. This relation may be useful in pointing systems 
which are evolved and/or reddened as the two effects move stars away from the Short-Period Blue 
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Envelope (SPBE) which has a meaning somewhat similar to the Zero- Age Main Sequence hne in 
the color-magnitude diagrams. It may also uncover inconsistencies and misclassifications, such as 
inclusion of pulsating stars. The approximating relation shown for the SPBE, (B — V)spbe = 
0.04 P~^-^^, with the period in days (Rucinski 1998b), does not have any physical significance, but 
is used for simplicity. We used the observed B — V color indices, without any corrections for the 
interstellar reddening, which should be small for most of the nearby objects of the sample. 

The B — V indices used in Figure 6 have been corrected for the third-light contamination only 
for 44 Boo where the brighter companion dominates all observed characteristics of the close binary; 
{B — V)b = 0.94 (Fabricius k, Makarov 2000). We were unable to evaluate corrections to the color 

indices for the five multiple systems which remained in the final sample. For VW Cep, KP Peg, 
V2388 Oph and ES Lib, the light contribution of the third body is moderately small {fl < 0.22), 
while for AA Get, with (5 = 0.59, we have no information on the color index of the companion, but 
its spectral type is similar to that of the binary (Tokovinin 1997). 

As can be seen in Figure 6, two objects have abnormally blue color indices, V445 Cep and 
KS Peg. With no spectroscopic data currently available for these stars, these peculiar colors have 
no explanation. We cannot exclude a possibility that these are in fact pulsating stars, although the 
well-defined double-wave light curve of KS Peg looks very typical for a close binary system. 

The absolute magnitudes derived from the HIP parallaxes are listed in Table 1. They were 
determined using: My = V'^^^ + 51og7r — 10, where the parallax tt is in milli-arcsec (mas). They 
can be compared with the predictions of the calibration established by Rucinski & Duerbeck (1997), 
Mv{cal) = —4.44 logP + 3.02 [B — V) + 0.12. This calibration was judged to provide predictions 
with the mean standard error a{Mv) — 0.22. The comparison is shown in Figure 7 with the 
My {HIP) errors estimated from the parallax errors, as listed in Table 1 and with the My(caZ) 
errors uniformly assumed at 0.24 as a combination of the calibration random error and the color 
index uncertainty. 

The My calibration is in perfect agreement with the data for the EW systems, but some 
deviations are clearly present for EB binaries and ELL variables. The deviations are expected for 
the EB systems as the color index of the dominant companion is not necessarily the best measure 
of the combined surface brightness of both components. In fact, the calibration is surprisingly well 
behaving even for these systems. The only case of a large deviation is the system AA Get which 
has a very poorly determined parallax, obviously due to a coupling with the visual-system orbital 
motion. Since the system with the observed spectral type and the short period cannot be as distant 
as luminous as implied by the measured parallax {My {HIP) = 0.41), the value of My(caZ) = 2.47 
has been adopted in this case for further use. This is the only system in the sample for which such 
a substitution was deemed necessary for a proper evaluation of the luminosity function. 

Two ELL systems show large and significant discrepancies between the two estimates of My, 
V335 Peg and IW Per, in both cases My {HIP) implies a fainter system. We have no explanation 
for these deviations. 
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8. THE LUMINOSITY FUNCTION 

The final sample consists of 32 objects. Counting in one-magnitude wide y™"^ bins centered on 
ymax _ 5^ numbers are: 2, 6 and 24. This progression with the apparent magnitude is very- 

close to the expected uniform spatial-density slope of 4x per magnitude and strongly suggests that 
our sample is indeed a complete one. We note that this rate of increase implies about 100 binaries 
in the next bin, 7.5 < y""^' < 8.5. This is not observed: the total number of EW+EB+ELL 
systems in the combined HIP and GCVS catalogs is currently 53. Thus, our choice of limiting the 
sample to V™'"'^ < 7.5 appears to be a correct one. 

The number statistics of the My, as listed in Table 1, can be used to derive the luminosity 
function by taking into account the changing volume for each increment of My. The luminosity 
function calculated in this way is given in Table 2 and is shown in Figure 8. The errors of the 
luminosity function have been estimated by scaling the Poisson errors by the appropriate volume. 
While the faintest bin, +4.5 < My < 5.5, contains two systems, the bin -1-3.5 < My < -1-4.5 
is empty. When plotting the luminosity function in Figure 8 in logarithmic units, we arbitrarily 
entered one system for this bin which is consistent with the implied Poisson error for such a 
substitution (1 ± 1 systems); this substitution was not used in the spatial density estimates below. 
Paucity of binaries in the luminosity interval which is entirely accessible to detection of the stellar 
variability strongly illustrates the strong influence of the Poisson fluctuations for faint systems of 
the sample. 

The combined spatial density for all EW-|-EB-|-ELL systems, p{all) = (1.02 ± 0.24) x 10~^ 
pc~^, is basically identical to the determination of Duerbeck (1984). The uncertainty of this 
estimate has been evaluated by adding the volume-scaled Poissonian error contributions from all 
My bins. The density estimate remains basically the same when very few bright systems in the 
range —0.5 < My < +1.5 are excluded, p(> +1.5) = (1.01 ± 0.31) x 10~^ pc~^, and is reduced 
only slightly when one considers only the fainter systems, p{> +2.5) = (0.90 ± 0.38) x 10~^ pc~^. 
Thus, the spatial density of contact binaries is dominated by faint systems, for which the Poisson 
statistical errors in the current sample are very large. The current uncertainty of the spatial density 
estimate is at the level of a factor of about two times. 

The Main Sequence luminosity function (Wielen et al. 1983) is shown in Figure 8, scaled by 
factor of 500. Such a scaled MS function is consistent with the new data for luminosity bins with 
My > +1.5. Table 3 gives the detailed comparison of the current sample with the scaled MS 
luminosity function. The predicted numbers of contact binaries for the RFO of 1/500 are listed in 
Table 3 under Np^ed- The next column in the same table gives the ratio of the actually observed 
numbers, No^g to Npred- The large Poisson errors of the predicted numbers indicate that the scaling 
by the factor of 500 is acceptable for My > +1.5, with the empty bin of +3.5 < My < +4.5 being 
a 2-sigma fluctuation. For the brighter systems within +0.5 < My < +1.5, a further scaling factor 
of about 10 would appear to be applicable, while for —0.5 < My < +0.5, the additional scaling 
factor would be about 60. 



- 12 - 



The assumed value for the RFO of 1/500 requires an explanation. In fact, the ratio of the 
luminosity functions, LF(E)/LF(MS), where E signifies all binaries considered here, gives the RFO 
between 1/400 and 1/900 (sec Tabic 3, the empty bin +3.5 < My < +4.5 is excluded). Only 
when an global average is calculated, one may recover Duerbeck's value of 1/1000 by inclusion 
of more luminous objects^. The spatial density of contact systems plunges down for My < +1.5 
(spectral types around A2 - A4) to the RFO ~ 1/5,000, dropping further to some 1/30,000 for 
My < +0.5 (spectral types around B8 AO). Although the number statistics is poor, the very 
low RFO for short-period binaries at high luminosities is a certainty. The very low spatial density 
of the bright contact systems was noticed before (Rucinski 1998a, b), but - in the contex of this 
paper - is actually expected given the assumed period limit at one day which naturally eliminates 
all massive and large binaries. 

In addition to the luminosity function derived for the current sample and the scaled MS 
function. Figure 8 shows the luminosity function derived previously from the 3 kpc deep Baade's 
Window sample of Disk Population contact systems using the OGLE-I data (Rucinski 1998b). 
We can clearly see a major discrepancy in the density estimates by factor of a few times over the 
whole range of luminosities. Since the RFO was estimated before at the level of 1/130 and now we 
estimate it at about 1 /500, the discrepancy involves a large factor of about four times. We suspect 
now that the OGLE-I sample either may have suffered from an observational problem or that the 
discrepancy tells us something very important about the difference between the local population 
and the one in the Galactic plane, at distances of 2 to 6 kpc. We return to this matter in Section 11. 

9. THE PERIOD FUNCTION 

The discussion of the orbital period distribution for the combined sample of open cluster and 
Baade's Window contact systems (Rucinski 1998b) suggested a strong concentration of such systems 
in the period range 0.25 to 0.5 days, with the peak around 0.35 days (see Figure 9). The peak 
was found to be even narrower and more concentrated at short periods for the binaries in globular 
clusters (Rucinski 2000). It may be argued that the short-period concentration is an artifact of the 
use of the linear units because - in this representation and for a flat log P distribution (see below) 
- the shorter-period intervals contain relatively more objects. However, the linear representation is 
ideal in demonstrating the abruptness of the well established cutoff at 0.22 days. We note that this 
cutoff still remains entirely unexplained, in spite of several efforts (Rucinski 1992b; Stcpich 1995; 
St§pieh et al. 2001). Most probably, the cutoff is related to the angular-momentum loss through 
magnetic activity, but - paradoxically - we do not know which of the two extreme assumptions is 
the correct one, i.e. whether the absence of late-type, short-period contact systems (Sp>K5, < 0.22 
days) results from the short period binaries being too active to exist for long time with short periods 



*If the RFO were really 1/500 for &4v > +1.5, the spatial density of the contact systems would be p(l/500) = 
1.51 X 10~^ pc~^, i.e. by 50 percent higher than the current estimate, but consistent with its error. 
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or not active enough to appreciably shorten their periods in the Hubble time. 

With the current sample, we cannot analyze the period distribution directly, by simply counting 
stars, because of the underlying correlation between the orbital period and luminosity which requires 
taking into account the varying search volume. In Rucinski (1998b), a novel concept of the "period 
function" (PF) was introduced. This function is basically the same as the period distribution except 
that the number of systems is expressed per unit of volume. The period function could be evaluated 
for the new data by summation of the contributions in "layers" of My, with the appropriately scaled 
Poissonian errors of the sums. The results for the current sample are given in Table 2, in both 
linear and logarithmic units of the period. The functions are shown graphically in Figures 9 and 
10. In both figures, the results are compared with those for the Baade's Window OGLE-I 3 kpc 
and 5 kpc sub-samples. We recall that the 3 kpc sample was better defined for the low-luminosity, 
short-period end of the period distribution, while the 5 kpc sample was preferable for P > 0.6 days 
because of the improved statistics for rare high-luminosity objects. It is striking how different the 
OGLE-I and the current samples are. Except for the shortest-period bin (0.2 < P < 0.3 days) 
and the longest period range (> 0.7 days), where the PF's agree within errors, the discrepancy is 
typically at the level of ten times. The missing systems in the +3.5 < My < +4.5 luminosity bin 
may be the main reason for the depressed PF, but the discrepancy is everywhere large and mainly 
reflects the very different estimates of the spatial density from both surveys. 

Except for the overall scaling relative the OGLE-I results, the period function derived here is 
interesting in that it does not show the short-period peak. Bearing in mind the poor statistics at 
short periods, the function is in fact consistent with an assumption of its flatness in logP units, 
which would indicate continuation and extension of the period distribution as seen at longer orbital 
periods. The seminal study of Dequennoy &; Mayor (1991) showed that the period distribution of 
all binary systems, over the very wide range of periods can be represented by a Gaussian cTirve in 
logP, with a shallow maximum at log P(d) ~ 4.8. As pointed by Heacox (1996), the Gaussian part 
introduces only a small curvature on distribution which is dominated by the logarithmic component, 
f{P)dP = P~^dP. Thus, the flat distribution /(logP) would appear to be appropriate and in 
fact seems to be valid for the contact systems discussed in this paper within the relatively short 
interval of 0.25 < P < 1.0 days. The main point is that we see no additional "piling up" at the 
very short-period end of this interval, as was suggested in Rucinski (1998b). 

10. THE AMPLITUDE DISTRIBUTION 

As was shown in Rucinski (2001), the distribution of the photometric variability amplitude 
is a useful tool for studies of contact binaries. In the present context, it relates to the detection 
constraints for small-amplitude systems; it is also dependent and carries information on the mass- 
ratio distribution of contact binaries. The examples used in Rucinski (2001) included the OGLE-I 
sample as well as a preliminary version of the sample discussed in this paper. The latter was based 
on objects selected on the basis of the original Hipparcos Catalog classification and in this respect 
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differs from the current sample which underwent a much more careful scrutiny for membership and 
variability type classification. 

The amplitude distribution (after correction for presence of third components in some systems) 
for the combined EW+EB+ELL sample is shown in Figure 11. The sample distribution is compared 
with the calculated amplitude distribution for the fill out parameter / = 0.25 and for a flat mass- 
ratio distribution, Q{q) = const. In contrast to the OGLE-I sample, which required a skewed 
Q{q), with a strong preference for small q, the current sample does not contradict the possibility 
of a flat mass-ratio distribution, although the statistics is obviously very poor. This discrepancy 
in the amplitude distributions can be taken as one of the strongest indicators that the OGLE-I 
sample was subject to image blending which systematically lowered the amplitudes and influenced 
the distance determinations. 



11. WHAT WENT WRONG WITH THE DENSITY ESTIMATE FROM THE 

BAADE'S WINDOW OGLE-I SAMPLE? 

The large discrepancy between the current estimate of the Relative Frequency of Occurrence 
of about 1/500 (for My > +1.5) compared with the RFO derived from the OGLE-I data (Rucinski 
1998b), as was found in Section 8, requires an explanation. While the variability-detection com- 
pleteness is a strong asset of the current sample, the weakness is in its small number statistics. 
As we suggested before, one direct and largest manifestation of this weakness is absence of any 
objects in the important luminosity bin of +3.5 < My < +4.5, which we suspect is a 2-sigma 
Poisson fluctuation. The current determination of the RFO is in an agreement with the estimate 
of Duerbeck (1984) who evaluated the average value at 1/1000. The results are consistent when 
allowance is made for the inclusion of intrinsically bright, but rare, systems in his sample. The 
close agreement in the average spatial density at about p = 1.0 x 10^^ pc~^ may be fortuitous, but 
- most likely - is a strong confirmation that proper corrections have been applied in Ducrbcck's 
determination for the low-inclination, low-amplitude systems which were only later discovered by 
the Hipparcos mission. 

Of concern is the large discrepancy of the present estimate of the RFO with the OGLE-I based 
result of RFO ~ 1/130. Here we list the several reasons why the Baade's Window sample has given 
a different - and very high - estimate of the RFO. 

1. For a pencil-beam survey, an error in the distance scale enters into a density estimate in 
the third power. The My (log P, B — V) calibration still requires refinements, especially for 
extreme colors indices and periods. The Malmquist bias corresponding to the standard error 
of the calibration, a{My) = 0.22, is also moderately large; for a uniform spatial distribution 
it is expected to lead to a systematic error in density of about 30 percent. 

2. While the current sample is strictly local, with the depth of some 40 to 400 pc, depending on 
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the luminosity range, the OGLE-I survey included all contact binaries of the Disk Population 

toward Baade's Window within a wide range of distances from about 2 to 8 kpc. The nearby 
systems were excluded by the bright limit of the OGLE-I photometry. All statistical inferences 
presented in Rucinski (1997a,b, 1998a,b) were based on two sub-samples to 3 kpc and 5 kpc 
in depth, but without any nearby systems which were cut off by the bright limit at around 2 
kpc. 

3. Blending of images in the crowded area of Baade's Window is expected to result in system- 
atic underestimates of distances. Of particular importance would be underestimation of the 
distances for systems located beyond the respective limits at 3 kpc and 5 kpc resulting in 
shifts into the respective sample volumes. The blending hypothesis is consistent with the ex- 
planation of the abnormal photometric amplitude distributions for the two OGLE-I samples, 
which are practically devoid of large-amplitude variables (Rucinski 2001). Blending of images 
is very hard to quantify without access to the original data. We note that even if all binaries 
were blended with equally bright stars (a shift in m — M = 0.75 mag.), the increase in the 
sample depth would be by 1.35 times, whereas the spatial-density discrepancy demands a 
factor of 4^/^ ~ 1.59. However, simple scaling may be entirely inappropriate as the numbers 
of contact binaries were seen to increase dramatically with the distance so that moderately 
small errors at the sample far limits may have played an enhanced role in density estimates. 

4. The Relative Frequency of Occurrence was utilized in the OGLE-I data interpretation. While 
it is a very useful concept, especially when numbers of variable stars can be directly related to 
numbers of stars subject to photometric monitoring (as in stellar clusters), the total number 
of monitored stars was not easily available in the particular case of the OGLE-I search, mostly 
because of the heavy crowding in this area. Evaluation of the RFO had to be done only in 
reference to the predicted density of Disk Population stars all the way to the Bulge (Rucinski 
1998b) , which in turn depends on the assumed value of the galactic exponential length scale, 
hR, currently a rather poorly known quantity. 

5. While the level of contamination by misclassified Population I pulsating stars should be the 
same for the current and the OGLE-I samples^, the one caused by the Population II variables 
of the RR Lyr-type and SX Phe-type may be larger in the direction of the Galactic Center. 
One can expect that the number of pulsating stars should be always smaller than the number 
of contact binaries (because the former exist only within the instability strips while the latter 
can be made from any MS stars), but this type of contamination cannot be a priori discounted. 
The local density of the RR Lyr stars is 6.2 x 10~^ pc~^ (Suntzeff et al. 1991), i.e. about 
1,600 times lower than the one derived here for contact binaries. An increase in the spatial 
density following r^^ or r^^ laws would produce an increase by an order of magnitude or so, 
hence the RR Lyr contamination is not expected to be very strong. 



^However, the current sample has an advantage of potential spectroscopic support in weeding out the pulsating 
variables. 



-16- 



6. Finally, we cannot exclude a possibility of a genuine increase in the number of contact binaries, 
especially the short-period ones, in the inner parts of the Galaxy. The globular cluster sample 
(Rucinski 2000) shows a period distribution even more strongly skewed to short periods than 
the Old Disk OGLE-I sample, with a narrow peak around 0.25 - 0.5 days. Apparently such 
systems are not very common in the solar neighborhood. The inner Disk may contain a 
richer admixture of very old, short-period, low-luminosity systems, similar to those found in 
globular clusters. 

None of the above causes is expected to lead to the systematic error in the space density by a factor 
of at least 4 times; however, it is quite possible that several of them actually combined in elevating 
the OGLE-I estimate. 

12. FUTURE WORK 

While the present, 7.5 magnitude-limit sample still suffers from several problems, it is a good 
starting point for establishing a volume-limited sample by extending it deeper in the magnitude 
scale. We noted that an improvement in the depth by one magnitude, to the limit of V^'^^ = 8.5 is 
expected to result in a substantial increase in numbers by about 100 objects. Such objects should 
be discovered through searches among bright stars advocated by Paczynski (1997, 2000, 2001). We 
have found that the current catalogs contain only about one half of the binaries which should be 
detectable within the one magnitude interval beyond the current limit, 7.5 < F"*"^ < 8.5. It is 
interesting and instructive to consider what will be these stars if we indeed increase the limiting 
magnitude to V"^""^ = 8.5. Simply scaling up by factor of four times, wc can predict that most 
system to be discovered will be moderately luminous, moderately long-period (0.4 to 0.7 days), late- 
A and early-F type binaries within +1.5 < My < +3.5. Some 60 - 70 binaries of those new 100 
systems will belong to this group, among them many well-known favorites, including the prototype, 
W UMa. The brighter and fainter systems will be more difficult to discover: At the bright end, 
for My < +1.5, an immense volume opens up very bright contact binaries, but they seem to be 
genuinely rare in space. We see a strong decrease of the spatial density in the present data, which is 
partly caused by the adopted period limit of one day eliminating massive, large systems, even when 
on the Zero Age Main Sequence. Rapid evolution of such systems away from the ZAMS additionally 
eliminates compact, early-type systems. The OGLE-I data suggest a particularly abrupt cut-olf 
for luminous systems with P > 1.3 — 1.5 days (Rucinski 1998a). Thus, we can expect about 10 - 
15 binaries with My < +1.5. At the other, faint end, with My > +3.5, probably also few systems 
will be detected. Although such systems are the most abundant in space, the search volume will 
remain small, so that deepening of the survey by one magnitude is expected to yield only some 15 
- 20 intrinsically faint systems. 

The implications of these considerations for the best strategy of establishing the spatial density 
of contact systems are as follows: (1) The pencil-beam searches are most promising, but they must 
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be free of the photometric blending; (2) Wide-sky surveys are not expected to give immediate 
results, if conducted to shallow limiting magnitudes. However, since periods and luminosities 
correlate for contact binaries, it may be possible to conduct targeted (with period constraints) 
deeper, wide-field searches. This may be the best strategy in the sense of a possibility of applying a 
combined, photometric and spectroscopic approach to learn most about these stars through detailed 
and full analysis. 

This paper is expected to provide a departure point for further studies. A spectroscopic 
survey to the same limiting magnitude would provide a uniform material to derive full physical 
parameters of the systems (via Ki radial- velocity semi-amplitudes). The mass-ratio parameter, 
q, which is crucial for light-synthesis solutions, would be then securely determined as the ratio of 
the Ki values. The radial-velocity orbit data would enable a study of spatial velocities, through 
combination of the Hipparcos tangential velocities with the center-of-mass velocities, Vq. Some 
efforts has been started in this direction at the David Dunlap Observatory where 60 radial velocity 
orbits, with the accuracy in Ki, K2 and Vq at the level of about 1 km s""-*^, have been recently 
determined (for complete references, see Rucinski (2002)). This spectroscopic survey very much 
needs an extension into the Southern skies. Availability of the mass-ratios will permit to address the 
currently missing mass-ratio term in the absolute-magnitude calibration. My = Mv{iogP, color, q) 
(Rucinski 1994; Rucinski & Duerbeck 1997). 

Spectroscopic observations are needed to check the variability classification of 14 ellipsoidal 
systems in the present sample. Currently only 6 of them have any spectroscopic data. There is a 
possibility that some ELL variables may turn out to be pulsating stars. Gerald Handler (private 
communication) estimates that about 30 percent of A/F type Main Sequence stars are S Set stars 
(the same percentage is expected for SX Phe Population II stars, but these are very rare in the solar 
neighborhood). The S Set stars have very short periods, when compared with contact binaries of the 
same type (Handler &: Shobbrook 2002), so that they are relatively easy to identify with combined 
spectroscopic and photometric data. However, variables belonging to the new class of 7 Dor (Kaye 
et al. 1999; Handler 1999; Henry et al. 2001; Handler & Shobbrook 2002), with pulsation periods 
from about 0.4 days and extending to 3 days, present a great danger of being taken for ELL binaries. 
Fortunately, most of them show several pulsation modes simultaneously excited. 

The current situation with spectroscopic data is summarized in Table 1. Systems marked 
SB2 are double-lined ones which currently have sufficient spectral information. Systems marked 
SBl, with one set of lines in the spectra, should be re-obscrvcd because - for short-period systems 
- it is very rare that modern observing techniques would not show a spectral signature of the 
fainter component. All systems marked by a colon should be also re-observed because the available 
radial-velocity data are inadequate. 
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IS. CONCLUSIONS 

The current sample of contact (EW) and related (EB and ELL) systems with orbital periods 
shorter than one day and selected from among bright stars to V"^"-^ = 7.5 (with a large frac- 
tion detected by the Hipparcos satellite) is expected to be complete in the sense that it contains 
all binaries showing photometric light variation larger than about 0.05 magnitude. This level of 
photometric accuracy implies that - according to predictions of Rucinski (2001) (Section 3.3 and 
Figure 4 there) - only about 10 - 15 percent of variables remain undiscovered due to low orbital 
inclination angle. Therefore, even if current statistics do not take these undetectable systems into 
account, the spatial-density estimate of p = (1.0 it 0.3) x 10~^ pc~^ is expected to be sound and 
consistent with the Relative Frequency of Occurrence, RFO 2± 1/500, when compared with the 
Main Sequence stars with My > +1.5. While both quantities, p and RFO, are still uncertain by 
a factor of about two times, mostly because of the small number statistics, the estimates should 
be relatively free of systematic errors. The spatial density estimate derived here is very close to 
that of Duerbeck (1984), but does not agree with the high value estimated on the basis of the Disk 
Population Baade's Window sample based on the OGLE-I data. As discussed in Section 11, while 
the discrepancy is most likely due to observational effects which produced the high value for the 
OGLE-I sample, we cannot exclude a possibility that the population of contact binaries observed 
towards the Center of the Galaxy is different from the one in the solar neighborhood. 

Special thanks are due to Hilmar Duerbeck, Gerald Handler, Hal McAlister and Bill Hartkopf 
for very useful suggestions and comments. 

Support from the Natural Sciences and Engineering Council of Canada is acknowledged with 
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A. NOTES ON INDIVIDUAL SYSTEMS 

The references cited below should be complete as of the end of the year 2001. 

44 Boo. A very well-known, frequently-studied system, the fainter component of the visual binary 
ADS 9494. Another common name i Boo (preferably the names should not be combined into 
44i Boo). The most general discussion in Hill et al. (1989a). Many speckle-interferometry studies. 
The newest, more accurate radial velocity orbit in Lu et al. (2001). The data on third-body light 
contribution from Fabricius &z Makarov (2000). The spectral type for the binary is estimated at 
K2V from the color index B — V = 0.94. The combined spectrum is dominated by the brighter 
GIV companion. 

VW Cep. A well-known, frequently-studied system. The most extensive study with references 
to earlier studies, Hendry & Mochnacki (2000). The data on third-body light contribution from 
Fabricius & Makarov (2000). 
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V759 Cen. Fainter at light maximum than 7.5 mag. hmit in iJp, but shghtly brighter in V, in 
accordance with the discovery data of Bond (1970), V'^"'^ ~ 7.45. Lacks a radial velocity study. 
The HIP light curve poor, possibly light curve changes. 

HT Vir. In a tight (0.56 arcsec) visual triple system, the slightly brighter component of ADS 9019. 
Excluded from the sample since an allowance for the third star has pushed its V"^"'^ to 7.86. A 
very interesting system, with a good radial- velocity orbit in Lu et al. (2001). 

KR Com. Discovered by Hipparcos. No reason for the type EB as in the HIP catalog. A radial 
velocity study in Rucinski et al. (2002). Many speckle- interferometry studies of the visual binary 
ADS 8863. Excluded after allowance for the light of the third star. 

V566 Oph. A well-known, frequently-studied system. An extensive spectroscopic study with 
earlier references in Hill et al. (1989b). 

VW Pic. Discovered by Hipparcos. A triple system with the fainter (1.9 mag.) visual companion 
separated by 9.9 arcsec. Called EB in HIP, but the eclipses are of equal depth. The Fourier 
decomposition of the light curve gives the type EA, so the this system has been excluded from the 
sample. 

AW UMa. One of the most frequently studied W UMa-type systems, mostly because of its 
extreme mass-ratio of q = 0.08 and presence of total eclipses. Recently relegated to the second 
place in the extremity of the mass-ratio, in favor of SX Crv (Rucinski et al. 2001) with q = 0.066. 
The most recent spectroscopic data in Rucinski (1992a). 

V972 Her. Discovered by Hipparcos. Called EB in HIP, but the light curve too shallow for 
classification, so classified as ELL. The spectroscopic orbit (Rucinski et al. 2002) indicates a low 
inclination systems. 

V445 Cep. Discovered by Hipparcos, but the period is in fact two times longer than given in HIP. 

CN Hyi. Discovered by Hipparcos. Included in Fabricius &; Makarov (2000), but only one visual 
component given, so most probably not a visual triple. A nice, well defined HIP light curve. 

MW Vir. Discovered by Hipparcos, but the period in fact two times longer than given in HIP. 

V867 Ara. Discovered by Hipparcos, but the HIP light curve rather poorly covered. A member 
of a visual binary (Fabricius &; Makarov 2000). Excluded from the sample because the Fourier 
decomposition of the light curve suggests the variability type EA. 

KS Peg. Called EB in HIP, but the light variation too small to classify the type. Hube & Gulliver 
(1985) and Hube et al. (1988) presented spectroscopic and photometric observations and suggested 
an extreme mass ratio, q < 0.1. The currently ongoing DDO radial velocity observations encounter 
difficulties as the broadening function is very wide, but does not split. The system requires an 
in-depth analysis. 

AA Cet. Discovered by Bloomer (1971) who classified the binary as EW and found the secondary 
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eclipse total. The brighter component of a visual binary ADS 1581 with the separation of 8.4 
arcsec (Tokovinin 1997). The (first published) light curve obtained by Hipparcos is definitely of an 
EB binary. This bright binary requires both photometric and spectroscopic study. The secondary 
eclipse may be total. 

V918 Her. Discovered by Hipparcos. No other data currently available. Called EB in HIP, but 

included among ELL here. 

V353 Peg. Discovered by Hipparcos. No other data currently available. Called EB in HIP, but 

included among ELL here. 

e CrA. The brightest y™"^ = 4.73 and one of the most frequently observed contact systems. The 
HIP light curve is rather sparsely covered. Light curve studies of Tapia (1969), Hernandez (1972), 
Twigg (1979) gave Qphot = 0.113 ± 0.002, while the spectroscopic study of Goecking & Duerbeck 
(1993) gave Qspec = 0.128 it 0.014, so no discrepancy in the mass ratio here, confirming good quality 
of data. 

RR Can. Major photometric studies: Knipe (1961), Knipe (1965), Chambliss (1971), Mochnacki 
&; Doughty (1972), Twigg (1979) gave consistent results, mostly because of the presence of total 
eclipses. The most recent spectroscopic study by King & Hilditch (1984) could be improved because 
of the low spectral resolution used. 

IS CMa. A large amplitude W UMa-type system discovered by HIP. The system begs further 
work. 

V851 Ara. Discovered by HIP. A shallow and sparse HIP light curve. No other data available 
except being included in the search for A Boo stars (Pauzen 2001) with negative result. Called EB 
in HIP, but the light variation too small to classify the type, so called ELL. 

LL Vel. The variability type uncertain (Matthews 1990; Krisciunas & Handler 1995; Mantegazza 
& Poretti 1995) because of the very small light variations, so ELL as in HIP is confirmed. 

V535 Ara. Photometric studies by Chambhss (1967) and Schoffel (1979). An UV light curve 
by Eaton (1991). The spectroscopic orbit of Schoffel (1979) certainly leaves a lot of room for 
improvement. 

S Ant. Not much work has been done on this nice EW system since the classic photometric and 
spectroscopic study of Popper (1956) where the binary was described as a single- lined system. It 
should be easy to determine a nice SB2 spectroscopic orbit with a modern equipment. 

VIOIO Oph. Photometric studies by Leung (1974) and Leung & Wilson (1977). Three SBl 
spectroscopic orbits summarized in Worek et al. (1988). An attempt to detect spectral lines of the 
secondary should be made. 

V449 Aur. Discovered by Hipparcos. Excluded as the HIP light curve suggests a detached binary 
with shallow eclipses. 
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2082 Cyg. Discovered by Hipparcos. No other data currently available. 

KP Peg. Discovered by Walker (1988) as a variable member of the visual binary ADS 14977 
(Tokovinin 1997). Light curve in Keskin & Can Akan (1989). 

PP Hya. Discovered by Hipparcos. No other data currently available. 

FO Vir. Discovered by Eggen (1983) and extensively observed in the following years (Schmidt 
& Fernie 1984; Antonello et al. 1984; Poretti et al. 1987). Mochnacki et al. (1986) estimated 
photometric mass ratio at 0.15, but the secondary has not been discovered spectroscopically yet. 
The spectral type of Gray & Garrison (1989) indicates that the star is slightly metal-poor. The 
star appears as CCDM J13298-I-0106A among visual binaries in the HIP database, but no other 
information is currently available. 

VI 130 Tau. Discovered by Hipparcos. The light curve looks like that of two strongly distorted, 
identical stars, but not of a contact binary, and this is confirmed by the Fourier decomposition 
of the light curve. The recent DDO spectroscopic data (Rucinski et al., private communication) 
indicate that it is a detached binary and definitely not an EB system. It is a weak-metal star (Abt 
1986; Gray & Garrison 1989). 

V2388 Oph. Rodriguez et al. (1998) discovered variability of the star, apparently independently 
of the HIP discovery. The star was a subject of many speckle interferometry investigations before 
discovery of the photometric variability. The revised visual binary orbit given by Soederhjelm 
(1999). The spectroscopic study of Rucinski et al. (2002) confirms the relative luminosity estimates 
of Soederhjelm (1999) and gives a good radial-velocity solution for both components. 

V335 Peg. Discovered by Hipparcos. No other data currently available. The photometric vari- 
ability very small so that the photometric period is uncertain and may require revision. 

TV Pic. The light curve shows a large O'Connell effect of the first maximum being higher then 
the second maximum. The photometric and high-resolution spectroscopic study of Pavlovski et al. 
(1998) discusses peculiarity of the system, but no radial-velocity orbit yet available. The binary is 
similar to TU Hor. 

ES Lib. A photometric study and a single-line radial-velocity orbit at low spectral resolution by 
Bartolini et al. (1973). The star has been a subject of several speckle interferometry studies. H. 
A. McAlister (private comm.) estimates that the visual companion must be faint, with Am > 2.5, 
which is confirmed by the apparent large ("undiluted") photometric amplitude of ES Lib. 

IW Per. Long recognized as an ellipsoidal variable (Morris 1985). Recent interest concentrated 
on the chemical peculiarity Adelman (1998); Pauzen & Maitzen (1998). Abt & Morrell (1995) list 
among the Am stars and estimate Vsini = 91 km s~^. 

DW Boo. Discovered by Hipparcos. No other data currently available. Glassified as EB in HIP, 
but this classification should be replaced by ELL in view of the very small light variations. 
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TU Hor. In the list of ELL variables of Morris (1985). The light curve shows a large O'Connell 

effect which very strongly contributes to the combined light variations. A photometric study by 
Ducrbeck (1977). The spectroscopic single-lined orbit by Duerbeck et al. (1979) [check] implies a 
low orbital inclination. 

GK Cep. Much studied system previously considered EW or EB, but excluded by us because the 
Fourier light curve decomposition indicates a detached binary. Pervious references in Niarchos et 
al. (1991). 

HL Dra. Discovered by Hipparcos. No other data currently available. The Fourier light curve 
decomposition indicates a detached binary, thus not confirming the EB type given in HIP. 

DX Aqr. The Fourier decomposition of the rather poorly covered HIP light curve indicates a 
detached binary. Member of the visual binary ADS 15562 Tokovinin (1997). A radial velocity 
study by Paffhousen &; Seggewiss (1976). 
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Captions to figures: 

Fig. 1. — Hipparcos light curves of eclipsing systems considered for the final sample {Hp magni- 
tudes). The orbital periods in days are given in parentheses by the names of the systems. The thin 
horizontal lines show the adopted maximum-light magnitudes H^""^ , which are given in numbers 
within the boxes. For those systems which did not have previously measured maximum-brightness 
V magnitude, these have been adopted for transformations into V^"'^; two systems fainter in Hp 
than 7.5 mag., V759 Oph and V566 Oph, have been retained in that process. Three systems, 
HT Vir, KR Com and VW Pic have been excluded because they were fainter than V"^"^' = 7.5 
when allowance was made for their companions. V867 Ara was found to be a detached binary, 
although its light curve is poorly covered and the EA classification may have to be revised. 

Fig. 2. — The same as in Figure 1 for the remaining systems. V449 Aur, V1130 Tau, GK Cep, 
HL Dra and DX Aqr were found to be a detached binaries. BB Scl and ER Vul, shown at the 
bottom of the figure, are the only two detached binaries fulfilling the period and brightness criteria 
with the available Hipparcos light curves. They are shown here for illustration purposes only. 

Fig. 3. — The Fourier coefficient a2 — 04 have been used to separate detached binaries (above the 
curve, triangles) from the EW (filled circles) and EB (open circles) binaries. The errors of the 
Fourier cosine coefficients have been determined by the bootstrap sampling process. In some cases, 
when the system happened to have a visual companion, the coefficients and their errors became 
magnified after allowance was made for additional light in the system. The open (EE) and closed 
(EW) symbols relate to the relative depth of eclipses, as defined in the next Figure 4. Some systems 
are labeled by their variable-star names. 

Fig. 4. — The Fourier coefficients 02 — ai used to separate EB binaries (open circles) from the rest 
of the systems. The other symbols are as in Figure 3. 

Fig. 5. — The Fourier coefficients 61 — ai used to find systems with light-curve maximum asymmetry 
(as measured by the first sine coefficient 61). The symbols used are the same as in the previous two 
figures. The observed correlation for ES Lib and VIOIO Oph suggests that these are semi-detached 
binaries transferring matter from the more- to the less-massive components, with an accretion spot 
on the side of smaller component. 

Fig. 6. — The period-color relation for the current sample. The symbols for the EW, EB and ELL 

systems are the same as in the previous figures. The continuous line gives the Short-Period Blue 
Envelope (SPBE) as derived in Rucinski (1998b); evolution and reddening arc expected to move 
stars down and to the right, away from the curve. Two ELL systems above the SPBE are V445 Cep 
and KS Peg. 

Fig. 7. — Comparison of the absolute magnitudes determined from the Hipparcos parallaxes with 
those estimated from the Rucinski & Duerbeck (1997) calibration. The symbols are the same as in 
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the previous figures. The error bars for My {HIP) are from the errors of the parallaxes while the 
errors in My(caZ) are assumed the same at 0.24 mag. 

Fig. 8. — The luminosity function for the whole sample (thick line) and for the sub-samples of ELL 
and EB systems (shaded and doubly shaded areas, respectively). The error bars show the Poisson 
errors in the bins. The Main Sequence function, scaled by factor 500, is shown by a dotted line. 
The OGLE-I estimate for the 3 kpc Baade's Window sample (Rucinski 1998b) is shown by the 
broken line. The error bars are based on Poisson statistics. Note that absence of low-luminosity 
EB and ELL systems suggested by this plot is not significant and may be explained by the low 
number statistics for faint systems. 

Fig. 9. — The period function for the combined sample (thick line) is compared with the data 
derived from the OGLE-I results (Rucinski 1998b). The broken line shows the period function for 
the 3 kpc sample (BW^), while the dotted line shows the same for the larger 5 kpc sample (SW5). 
The data for the latter are valid only for P > 0.5 days, but are statistically better defined than 
for the 3 kpc sample in this period range. The insert shows the long-period portion magnified by 
factor 5x (right vertical axis) to illustrate that, for long periods, the current period function is 
basically consistent with that for the 5 kpc OGLE-I sample. Note that linear units of the period 
are used in this plot. The next figure shows the same data in the logarithmic units. 

Fig. 10. — The period function for the whole sample (thick line), in the logarithmic units of 
the period. The sub-samples of ELL and EB systems are marked by shaded and doubly shaded 
histograms, as in Figure 8 (see the caption to this figure). The OGLE-I 3 kpc {BW3) and 5 kpc 
{BW5) samples (Rucinski 1998b) are shown by the broken and dotted lines. See the text and the 
previous figure for further explanations. 

Fig. 11. — The amplitude distribution for the sample. For systems in triple systems, the corrected 
amplitudes are shown. The ELL binaries occupy, basically by definition, the first bin at < 0.1 mag. 

because the criterion for ELL assignment was on the Fourier term, —0.05 < 02 < 0. The dotted 
line gives the calculated distribution for contact binary systems (Rucinski 2001) with the degree 
of contact / = 0.25 and the flat mass-ratio distribution, normalized to the same total number of 
objects. The shaded area gives the distribution for the OGLE Baade's Window sample to 5 kpc 
(which includes the 3 kpc sample), truncated to and normalized for objects with amplitudes larger 
than 0.2 mag. 
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Table 1. The 7.5 magnitude sample of EW+EB+ELL systems 
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Table 1 Contimiod 



Name 


HIP 


HD 


P 


H? 


EX? 


Type 


V' 


max 


B-V 


/3 


Sep 


ampl 


TT 


My 


10 


ai 


0.2 


a4 


bi 


RV 


Sp 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 
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+0.002 
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-0.001 
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H 


EX 


EA 


7.28 




0.19 






0.31 


6.4 ± 0.6 
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+0.005 
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-0.012 


+0.003 


SB2 


GOV 



Note. — Columns: (1) Name; (2) (3) HIP & HD; (4) period in days; (5) H = discovered by Hipparcos; (6) EX = excluded from the sample; (8) (9) V"^°'^ as observed and corrected for the third 
light; (10) observed B — V, uncorrected for the third light (see text re. 44 Boo); (11) /3 = L-i/L"}^^; (12) separation in arcscc; (13) observed amplitude of light variations in magnitudes, uncorrected 
for the third light; (14) Hipparcos parallax in milliarcsec; (15) My derived from the y""^ and the parallax; (16) - (19) the Fourier cosine coefficients, in light units, uncorrected for third light; (20) 
the sine Fourier coefficient; (21) current status of radial velocity observations; (22) estimate of the spectral type. 
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Table 2. Luminosity and period functions 



Mv 


Nobs 


dist 


LF 








N (logP) 












(pc) 


(pc-^) 


-0.65 


-0.55 


-0.45 


-0.35 


-0.25 


-0.15 


-0.05 


Luminosity function (LF) 


















-0.5 +0.5 




398.1 


3.78e-9 ± 3.78e-9 

















1 





+0.5 +1.5 


5 


251.2 


7.53e-8 ± 3.37e-8 














2 





3 


+1.5 +2.5 


18^ 


158.5 


10.08e-7 ± 2.54e-7 











2 


8 


5 


3 


+2.5 +3.5 


6 


100.0 


14.32e-7 ± 5.85e-7 








1 


4 








1 


+3.5 +4.5 


(1)' 


63.1 


(9.50e-7 ± 9.50e-7) 























+4.5 +5.5 


2 


39.8 


7.57e-6 ± 5.35e-6 





2 


















Period function (PF) 



total 
PF (e-7) 
crPF (e-7) 



2 
75.7 
53.5 



1 

2.4 
2.4 



6 
10.7 
4.8 



10 
5.1 
1.7 



6 

3.0 

1.3 



Period function (PF) - linear 



N (P) 



PF (e-7) 
aPF (e-7) 



0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


0.95 


75.7 


2.4 


10.7 


2.1 


4.2 


1.8 


3.7 


0.9 


53.5 


2.4 


4.8 


1.1 


1.6 


1.0 


2.5 


0.6 



Note. — The right head side of the LF section gives the number of systems in respective My and distance intervals 
in bins AlogP = 0.1, centered at logP (in days) as given in the headings. The PF section gives the resulting period 
function in the same intervals of AlogP(d) = 0.1 as well as a separately binned distribution of linear intervals of 
AP = 0.1 days. 

''Because of the large error of the parllax, AA Get has been moved from Mv(HIP) = 0.41 ± 1.11 to Mv{cal) = 
2.47 + 0.24. 

''As explained in the text, the interval +3.5 < My < +4.5 does not contain any systems so that the LF is given 
for the hypothetical case of 1 ± 1 system there. This assumption was not made when the total spatial density was 
calculated, as given in the text. 
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Table 3. Predicted numbers of EW+EB+ELL systems (RFO = 1/500) 



Mv 




Nobs 


pred 


Nobs/Npred 


LF(E)/LF(MS) 


RFO 


-0.5 + 


0.5 


1 + 1.00 


63.1 + 7.9 


0.016 


3.17e-5 


31,546 


+0.5 + 


1.5 


5 + 2.24 


55.5 + 7.4 


0.090 


1.80e-4 


5,291 


+ 1.5 + 


2.5 


18 + 4.24 


23.9 + 4.9 


0.75 


1.51e-3 


662 


+2.5 + 


3.5 


6 + 2.45 


10.8 + 3.3 


0.56 


1.12e-3 


892 


+3.5 + 


4.5 





4.90 + 2.2 








oo 


+4.5 + 


5.5 


2 + 1.41 


1.70 + 1.3 


1.17 


2.35e-3 


425 



Note. — The first four columns give a comparison of the number of observed 
systems, Nobs, with the predicted number assuming the Relative Frequency of 
Occurrence of 1/500, compared with Main Sequence stars of the same luminosity. 
The last two columns give the ratio of the luminosity functions and its inverse 
which is the RFO evaluated in individual bins. See Table 2 for a comment re. 
A A Get moved between the My bins. 



